Hall resistance in the hopping regime, a "Hall Insulator"? 
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The Hall conductivity and resistivity of strongly localized electrons at low temperatures and at 
small magnetic fields are obtained. It is found that the results depend on whether the conductivity 
. or the resistivity tensors are averaged to obtain the macroscopic Hall resistivity. In the second case 

Q\ ' the Hall resistivity always diverges exponentially as the temperature tends to zero. But when the 

Q\ , Hall resistivity is derived from the averaged conductivity, the resulting temperature dependence is 

sensitive to the disorder configuration. Then the Hall resistivity may approach a constant value 
5^ ' as T — > 0. This is the Hall insulating behavior. It is argued that for strictly dc conditions, the 

Oh, transport quantity that should be averaged is the resistivity. 
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Recently it has been stated that the zero-temperature Hall resistivity p xy (ui) of noninteracting electrons in 

the insulating regime remains finite as the frequency to — > 0. This result was derived from the Kubo formula for the 
frequency-dependent conductivity. It was found that at zero temperature the disorder-averaged a xy (u>) of the Anderson 
insulator vanishes at low frequencies proportionally to lo 2 . Since to leading order in w the longitudinal conductivity 
a xx ~ iujso, where e is the dielectric constant, this yielded that the Hall resistivity 

Pxy ^ ®xy /( cr L + cr L) approaches 

■ a constant in the small-frequency, zero-temperature limit. 

Obviously one would like to examine the Hall resistance at finite temperatures. Here we address this question for 
strongly localized electrons. We consider the problem using the Holstein model [Q for the Hall effect in a system 
with localized states. The conductivity tensor for a triangle of three sites is obtained first. To get the macroscopic 
Hall resistivity one may attempt to apply two averaging procedures. We find that the result strongly depends on 
whether the conductivity tensor is averaged before inverting it or vice versa. In the first case a Hall insulator behavior 
is possible, while in the second the macroscopic Hall resistivity increases with decreasing temperature. The latter is 

■ similar to the finding of Friedman and Pollak ||. The two behaviors are argued to be related with "ac" and "dc" 
| conditions. 

The previous discussions of the zero-temperature Hall resistivity were all based on the ensemble averaged 

' Kubo formula, which yields that a xx — a yy and that a xy vanishes proportionally to the magnetic field (the bar above 
. the quantity indicates ensemble averaging). However, before averaging, a xy includes a field-independent term. We 
show that in the strongly localized regime this term is comparable in magnitude to a xx and a yy . This leads to delicate 
cancellations when the local (unaveraged) conductivity tensor is inverted to obtain the resistivity, and consequently 

■ to the very different temperature dependences described above. 

It is convenient to investigate the transport properties of electrons in the hopping regime by constructing the rate 
equations for the electron distributions, utilizing the electronic transition probabilities between localized states. One 
has 

3 

in which rii is the electronic population of site i (the term 'site' is used for a localized state) and is the rate 
of the population decay by phonon-assisted hopping into site j. A delicate point is the dependence of the rate on 
the magnetic field, H. As was shown by Holstein Q, this is due to interference of the "direct" amplitude to hop 
from i to j, with the indirect amplitude via a third site £, i — > £ — > j. The magnetic field dependence of P+j then 
necessitates the consideration of at least three sites. To write this rate explicitly, we employ the Holstein model [|| for 
the electron-phonon interaction, in which the ion displacements are coupled to the local (site) density of the electrons, 
and denote by ti the single-particle energies of the localized states, which are assumed to be randomly distributed, 
and _by_ the overlap of two wave functions localized at sites i and j. The strong localization regime is characterized 



by [ffi 



Jij |<| etj = €i - Cj. One finds 
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Here P^ r arises from the direct amplitude alone and is independent of the magnetic field, 

/oo 
dt e" 12 ' e 2 ^(*)-9(0)) ) 
- OO 

q ^1 

where uo q is the phonon frequency, v q is the interaction matrix element and N q — l/(e l3uJq — 1) (% = 1). p™ dlr comes 
from the interference between the direct (1 — ► 2) and the indirect (1 — ► 3 — » 2) amplitudes 

i$* p = m(l - n 2 ) I Ji 2 J 23 J 3 i I 29 e* 32 

dtl e »(ti)+fl(t a )+»(ti+t a )-3fl(0) 

-OO 

(1 - n 3 )e" 12tl+ " 13t2 - n 3 e le i2*i+^32*2j ( ( 4 ) 
and contains the occupation of site 3, n 3 , and the magnetic phase, <^i 3 2, 

0132 = ~H ■ S, 
c 

S= (Ri x R3+R3 x R2 + R2 x i?i)/2. (5) 

Here -R, is the radius vector of site i and 5 is the vectorial area of the triangle. The field-dependent part of P12 
includes a term even in the field (proportional to cos0i 3 2) and a term odd in it (proportional to sin (^132). The first 
gives a correction to the direct rate and will be discarded henceforth. Obviously it is the term odd in the field which 
is capable of producing the Hall resistance. 

We now apply the rate equation ([j]) to a group of three sites, 1, 2 and 3, to obtain the current driven by an external 
ac electric field E of small frequency uj. In the presence of an electric field the occupations will be modified in a 
way that can be expressed by changes 8 fa in the local chemical potentials Q 

n i = n?-(3n?(l-n° i )5fa, (6) 

where nf is the Fermi distribution. Also, the rate Pij which depends on e.j — ej is changed to depend on e.j — e» + eE ■ 
(R — Rj). This way one obtains the response of the system to the local electrochemical potential differences |7j 

Cij = — — - — — ^ E ■ Rij, Rij = Ri — Rj. (7) 
Explicit calculations of Eqs. @ and (|J) yield 

e ^dT = e(P21 ~ Pl2 + Pil ~ Pia) 

= G 12 Ci2 + G 31 Ci 3 + sin0rC 23 , (8) 

with analogous expressions for the other two sites, where <f> is the magnetic flux enclosed in the triangle. These 
equations determine the electrochemical potential differences, Qj. However, in practice one does not have to solve for 
the C^s since Eq. (||) gives the current (see Eq. ( |l3| ) below.) In (||), Gy is the conductance of the bond ij, arising 
from the direct rate, 

Gij = Gji = 

/OO 
dte ieiit e 2(ff(t)-ff(0)). (9) 
-00 

At low temperatures and for weak electron-phonon coupling, the bond conductance becomes 

= e 2 /3 | Jij | 2 ^^ e-Kl^l+I^W^D, (10) 
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where AT is the phonon density of states. (Energies are measured from the Fermi energy.) The interference process 
leads to the phase-dependent term, with 

T = e 2 l3n\{\-nl) | J 12 J 23 J 3 i 

/CO 
dt 2 e 9( - tl)+9(t2)+9( - tl+t2 ^ 3g(0) 
-CO 

£l2tl ((1 - n°)e" 13t2 + n°e le32t2 y (11) 



Expanding this expression for weak electron-phonon interaction and using ( |10[ ) we find that T can be written in the 
form 

p _ | J12J23J31 I ( G31G12 



2^23^31 I / ^31^12 

4^0 V| J 31 J 12 |2n?(l-n°) 
G12G23 G23G31 



\ J12J23 \ 2 n° 2 (l - n° 2 ) I J23J31 | 2 n°(l-ng) 



(12) 



This expresses the fact that the indirect rate involves two scattering events by the phonons |0,p|. The temperature 
dependence of T, at low temperatures, can be obtained from (|l0|), using — ~ exp[— f3 \ |]. 
The current density is given by 

J=e£f*. (13) 

i 

If one now introduces the effective field E e tf which produces the electrochemical potential, 



Cij - -Eeff ■ R 
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E ef f = ^r[C2 3 i?i + C3ii?2 + CuRa] x z, (14) 



the current density becomes 



3 = vEeff 

a = R12R12G12 + R23R23G23 + R31R31G31 

+ T Sin ^23^1+^31^2+^12^3)- (15) 

A remarkable observation is that the part of j which is proportional to the magnetic field is perpendicular to the 
effective electric field. In deriving the expression for cr it was assumed that the triangle lies in the x — y plane, 
perpendicular to z. Note that E e ff, Eq. (|l4|), is invariant to the choice of the coordinate origin. 

The current response to the electrochemical potential difference is our central result. From (15) one finds for the 
conductivity tensor 

a xx = (R X 12 ) 2 G 12 + (R X 23 ) 2 G 2 3 + (i?3i) 2 G 3 i, (16) 
(pyy is given upon replacing Rfj by Rfj), and 

a f,l = R12R12G12 + -#23-^23^23 + ^31^31^31 



yx 



±2S z Tshx4>. (17) 



It is straightforward to check that these results follow very simply also from the Kubo formulation, with the 
understanding that it yields the response to the effective field, E e f f : 

ie 2 ui s—^ r< m\xi\n >< n\xj\m > 



< m\xj\n >< n\xi\m > 
UJ + e nm + ir\ 



(fm ~ /»), (18) 
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where 77 is a positive infinitesimal, |m > and \n > are eigenstates and f n , f m their populations. We generalize the 
derivation of j| to include electron-phonon coupling and consider the oj > 0, oj — > limit. For example, to get the 
leading, (l,2)-type terms, in (16) and in the first line of (17), the two relevant states are: |1 >: electron in site 1 
with a phonon bath in equilibrium and \2,q >: electron in site 2, with the same minus one phonon in state q, where 
oJ q = £21- Thus, the approximate eigenstates are: 



|m>=|l>+2^^^|2,<Z>, (19) 

V £l2 



|n>=|2,g>+^^|l>, (20) 
£21 — w„ 



where v q 
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< n\x\m >^ ±Z| q ' L \ (21) 



Using Eq. jl^), this produces the first term in jl^): 
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irfie Jf 2 (R12 ) ~2~ N q N (e 2 i ) \ u „ . (22) 

Using the expression similar to ( |2p| ) for the matrix elements of y, yields the first term in (|l7|). The (2,3) and (1,3) 
terms in Eqs. ( |l6| ) and ( |l7| ) are similarly obtained. To get the T term in (17), one has to mix in Eq. ( pj| ) also the 
state 1 3, q, q' > (Te. |3 > minus the phonons q and g') in two ways: 

(a) A straightforward correction Y^ g > [J\z v q v q' I (^13 + u q + u) q >)]\3, q,q' >, which is first order in J but second order 
in the electron-phonon interaction. It will be dominated by the "resonant" contribution (q' such that u q +ui q > — £31) 
and by the nonresonant contribution given by q' = q" such that w g » = £32- 

(b) The Holstein contribution, the mixing of | 2,q > via the intermediate state | 3,q,q' >: ^2 ql \JnJ-ii v q I v q' | 2 
/ (£13 + Lo q + u>q')(ei2 + oJ q )] 2, q > . Here, as found by Holstein, the resonant contribution with £31 = uj q + Lu q > will 
yield the needed phase to have a term odd in the magnetic field (the V term) in ([l7]). Putting all the above together 
we find within the required accuracy 



(n|a;|m} = 12 q R\ 2 — 2iTr5{e^i — to q — uj q i) 

UJ 

ri-i + rL-2 



J\zJz2 v q 



Rt 



(23) 



Using this in Eq. ( J18| ) produces the additional T term in Eq. (J17J). 

The result of the above calculations is the conductivity tensor a of a single triangle. We would like to have 
macroscopic quantities like the Hall resistivity p xy . For that purpose one might average tr over the orientations and 
sizes of the triangles and also over the on-site energies. A priori two averaging procedures exist. One can average 
(j, resulting in a, and then calculate p xy = (a~ 1 ) xy . One can also average the resistivity tensor of a triangle 
resulting in cr -1 and then calculate p xy = (cr -1 )^. As we will see these two procedures lead to qualitatively different 
results. 

We note that <j xy includes a term independent of the magnetic field, which is of the same order of magnitude as 
o~ xx - Were we to average a over directions before inverting it, then this term would have disappeared. However, if *p 
is to be averaged, then this term remains, and leads to delicate cancellations in the denominator of p. This, in turn, 
is the cause of the two different temperature dependences of the Hall resistivity. 

We first average the conductivity tensor over directions and then invert it. In that case, 

(1) 2S z rsin^ 

Pxy (WGV [ ] 

where we have assumed for simplicity that G12 is the largest conductance. To obtain the temperature dependence 
we consider the situation in which the magnetic field-free hopping conduction takes place along the bond 12 and site 
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3 supplies the interference path. Thus we imagine ei and e 2 to be below and above the Fermi level, but close to it, 
while €3 is away from the Fermi energy. Then [cf. Eqs. ( |Io|) and (|l2|)l 

P# ~ i ~ expH^a - 2e a + Cl )], (25) 

and is sensitive to the averaging procedure over the on-site energy distribution. One may imagine that the energy e 2 
is mostly in-between the energies t\ and e 3 , in which case the log-average of p^y will lead to a constant value for the 
Hall resistivity at very low temperatures - i.e. a "Hall-insulating behavior" . 

We next consider the transverse resistivity obtained by inverting the full conductivity tensor and then averaging 
over directions to restore rotational invariance. The result is 

ri\ 1 rsin</> 

p( 2 ) = _ , (26) 

xy 2S Z G12G23 + G23G31 + G 3 iGi2 



Due to the cancellations occuring when <r of Eqs. ( |16| ) and (17) is inverted, the denominator here includes G12G23, 

etc., but not the (larger) term G\ 2l a s in (p4|). Consequently, p£y increases exponentially as the temperature tends to 
zero, independent of the specific configuration of the single-particle energies. This is because of the factors n°(l — rq) 
in Eq. (n_2J) . Consider for example the energy configuration specified above. In that case the leading term in V is of 
order exp[— (3(e 3 — ti)] while G12G23 (G12 ~ exp[— f3(e 2 — ei)], G23 ~ exp[— fle^}) dominates the denominator in Eq. 
@), leading to p'iy ~ exp[/3e 2 ]. 

Both p^y and p^j are independent of the strength of the coupling to the phonons. This is in analogy with the 
"classical" (Boltzmann equation) result for the Hall coefficient, which does not contain the mean free path. 

The physically correct way of averaging may depend on whether the experiment is a dc one or an ac one. We 
consider a measurement to be an ac one if the frequency is such that the electron cannot traverse the sample from 
one current contact to the other during one field period. In this case the current contacts are irrelevant, the current 
is inside the sample, and the macroscopic current density is obtained by summing the contributions from all triangles 
within a unit volume. This summation is equivalent to the averaging of the conductivity tensor of a single triangle. 

We define a measurement as a dc one when the electron traverses the sample during a time short compared to the 
period of the field. The current is therefore flowing from one current contact to the other through a percolation chain 
of bonds. This means that the direction of the current in each bond is defined. To find the elementary Hall voltage 
created at this bond we have to invert the conductivity tensor of the triangle attached to this bond and express this 
elementary Hall voltage in terms of the resistivity tensor of a single triangle. The total Hall voltage is obtained in 
this case by summing over the bonds along the percolation chain. This summation is equivalent to the averaging of 
the resistivity tensor of a single triangle. 

To summarize: Two independent derivations of a at zero frequency but finite temperatures, were given for the 
Holstcin model. Here the Hall conductivity too has a finite dc value when real, phonon-mediated, transitions are 
allowed. Surprisingly enough, we find that the answer depends on which transport quantity is averaged over directions. 
It was shown that the ensemble averaging needed to get the macroscopic Hall resistivity is subtle and the result 
depends on whether cr is averaged before or after having been inverted. The former procedure leads to a possible 
"Hall insulating" behavior. The latter leads to a p xy which grows exponentially when the temperature is lowered. 
Using the percolating path picture of , one might speculate that the latter is the proper averaging procedure for 
the dc limit. 
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